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Tt ■ 1 Introduction 

The field of integrable systems is a vast one that has grown almost explosively 
in recent years. It is not remotely possible to review all the developments 
in the field in one hour, so I am restricting my attention to some particular 
ly-N ' topics. I shall only consider one dimensional quantum systems, or 1+1 as it is 

f^ , usually said, including the time dimension, but necessarily therefore including 

some discussion of the statistical mechanics of two dimensional lattice systems 
(vertex models) which share a common mathematical structure with the 1+1 
quantum systems, and have been the key to much of the understanding of the 
quantum systems. Q g [Q 

C^ • This is still far too big a field to review in any detail in one hour, so I shall 

G , focus on some generalizations of the spin one-half chain with nearest-neighbor 

coupling, for which Bethe first set down his celebrated ansatzQ. The simplest 
generalization, to fully anisotropic coupling, the so-called XYZ model, has an 
amazingly rich mathematical structure, parts of which are still being elucidated. 

^ , A second generalization, which has been very fashionable lately, is to extend the 

coupling beyond nearest neighbors, specifically to a long range inverse-square 
type coupling between spins one-half. This also turns out to be an integrable 
^^ ' system, although the wave functions are not of the Bethe ansatz type.[H Q 

'H How relevant are these models to real physical systems? The simplest 

isotropic spin one half XXX Heisenberg antiferromagnet is a good Hamiltonian 
for the quasi-one-dimensional system CPC, and one of the early experimental 
vindications of the Bethe ansatz was the confirmation by neutron scattcringH 
that it correctly predicted the observed elementary excitation spectrum, in con- 
trast to the standard spin wave theory used at the time. An appropriate con- 
tinuum limit of the XYZ chain gives the sine-Gordon model, and the XX Z 
chains are good representations of known systems. The closely related one- 
dimensional Hubbard model has been widely used to describe one-dimensional 
conductors, and may be relevant to some high temperature superconductors. 
One of the great successes of the Bethe ansatz has been the solution of the 
Kondo problemp. The Luttinger liquid model[pl, first used to analyze low 
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energy properties of the XXZ spin chain, is proving valuable in studying one- 
dimensional mesoscopic conducting sy stems |10[. In fact, Bethe ansatz models 
are very useful for one-dimensional conducting ring systems having a thread- 
ing magnetic field, because the Bohm-Aharonov phase can be incorporated as 
a twisted boundary condition which, it turns out, does not spoil the integra- 
bility. Another very fruitful interaction in recent years has been that between 
Bethe ansatz methods and conformal field theory, discussed below. Finally, the 
inverse-square systems [^ are cousins of the Bethe ansatz systems. They do 
not have Bethe ansatz wave functions, but their energy levels are given by Bethe 
ansatz like equations of a rather simple kind. These systems are closely related 
to the important edge states in the quantum Hall effect, and also, surprisingly, 
to level distributions in quantum chaotic systems |1 4 . 

2 Quick review of Bethe ansatz 

The simplest Bethe ansatz system to visualize is the "billiard balls on a line" 
problem — classically, imagine N perfectly elastic billiard balls on a ring, say, to 
give periodic boundary conditions. The initial set of momenta ki, ..., k^ is con- 
served in subsequent collisions. The quantum mechanical version, solved long 
ago by Lieb and Liniger|l3|, has N bosons interacting through repulsive delta 
function potentials. The two particle problem is equivalent to a single particle 
interacting with a static delta- function potential, and the standard boundary 
condition at the potential leads to a two-particle wavefunction which is a sum 
of two plane wave terms with the momenta permuted, 

i(;{xi,X2) = £^(''1^1+'=^^=) + e*^i2e^(''^^i+'=i"=^) (1) 

for xi < X2, and boson symmetry determines the wavefunction for xi > X2- 
The phase shift term is found to be: 
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iih - k2) ,2) 



c + i{ki - k2) 



where c is the strength of the repulsive delta-function potential. This is a simple 
Bethe ansatz wavefunction. Notice that it vanishes if ki — k2. 

Surprisingly, this wavefunction generalizes in a simple way to N particles, 

tp{xi,X2, ■■■,xn) = ^a(F)e'^''^3'^' (3) 

for xi < a;2 < ... < xn, with other orderings given by the boson symmetry. 
Here a{P) is a product of phase shift factors for each pair of transitions needed 
to get the ordering P from the ordering 1, ..., A'^. It follows immediately from 
the form of the phase shift term above that all the ki must be distinct, or 
the wavefunction vanishes identically. The total energy and total momentum 



are easily seen to be given hy E ~ ^ kf , K ~ ^ ki, from applying the E, K 
operators in regions where no two XiS coincide. 

The actual values of the ki appearing in the wavefunction are determined 
by the boundary conditions, just as they are for free particles. For periodic 
boundary conditions, the total change in the phase of the wave function on 
taking a particle around the system must be 27r/, where / is an integer. This 
phase is part kinetic, from the e^^^ term, and part phase shifting from passing 
the other particles so the total phase shift 

kJL + Y.^,^ = '^7^IJ (4) 

Now writing the phase shift in terms of the momenta, we get a set of coupled 
nonlinear equations for the fci's for a given set of quantum numbers li. Since 
the ki must all be distinct, the ground state is fermi-sea like, with a set of Ij 
equal to 0, ±1, ±2, ... . Elementary excitations are given by gaps (holes) in the 
sequence of integers, or single integers above the sea (particles). It is found 
that for an attractive delta function potential, the ground state for A'' particles 
is a string of pure imaginary momenta having uniform spacing ic. Actually, 
this attractive boson gas is a pathological system in the thermodynamic limit 
where all the particles bind together infinitely tightly. However, these strings 
of bound momenta prove to be a common (usually nonpathological) feature of 
Bethc ansatz systems, as discussed below. 

Let us now turn to the spin one-half XX Z chain having Hamiltonian: 

N 

where we shall restrict our attention to the gapless regime, |A| < 1, and write 
A ~ cos 7. The eigenstates of the Hamiltonian are described by Bethe ansatz 
wavefunctions having the same general structure as in the boson model, the 
bosons being replaced by magnons. A magnon is created by applying the opera- 
tor a^ = ^e*'^"5'+ to the all-spins-down "ferromagnetic" state. It is convenient 
to parameterize the magnon momenta kj in terms of so-called rapidity variables 

Vj-. 

g^fc, ^ sinh(t;j- + i-f/2) 

sinh(wj - ij/2) ^ ' 

The reason for the change of variables is that the magnon- magnon phase shift is 
a function of rapidity difference. It was not a function of momentum difference — 
the lattice spoils Galilean invariance. The equations are much easier to handle 
in the new variables. The boundary condition for allowed rapidities analogous 
to (4) (or, rather, its exponential) is: 

' iwMvjJrjnJTfX _ -A- sinh(t;j- -Vi + i-j) 
smh{vj — i^/2) J \^ sinh(wj — f ^ — 27) 



One significant difference from the (repulsive) boson system is tliat the 
magnon rapidities given by solving these equations include in general some 
grouped in strings of certain lengths, the allowed lengths depending on 7, and 
also showing a very complicated system size dependence for finite systems. 

One of the most important applications of the Bethe ansatz method is finding 
thermodynamic properties. Roughly speaking, this is done by summing over an 
ensemble of systems having the same macroscopic quantum number density with 
appropriate entropic weighting factors. For the spin chains, the presence of so 
many string-like excitations makes this a formidable task even in the infinite 
limit |14| , and for finite systems, where string lengths vary with system size, it 
is even more tricky. Fortunately, new techniques have recently been discovered 
which circumvent these difficulties [[la f|16|. We shall discuss them below. 



3 Luttinger liquids and conformal fields 

One real problem with the above analysis is that, although all the energy levels 
of the Bethe ansatz system can be found and the wavefunctions written down, 
it is very difficult in practice to find matrix elements of operators between these 
states, and hence connect the known excitation spectrum with what is observed 
in, say, neutron scattering from a spin chain. The first significant step in this 
direction was made by Luther and Peschelfll, who argued that the low en- 
ergy excitations of the gapless XX Z spin chain could be described using the 
Luttinger model, solved by Mattis and LiebpSJ. These ideas were extended 
and put on a firm footing by HaldaneH, who coined the term Luttinger liq- 
uid. Essentially, the method is to transform the magnons to spinless fermions 
using the Jordan- Wigner transformation, whereupon the XX term is seen to 
be a kinetic fermion hopping contribution to the Hamiltonian, and the A term 
is a four-fermion interaction. The fermion energy-momentum curve is then 
linearized near the "fermi surface" so that density fiuctuations can be repre- 
sented by bosons. With this approximation, the whole low-energy Hamiltonian 
is quadratic in the boson operators, and can be diagonalized and solved. The 
fermion creation operator is found to be closely related to an exponential func- 
tion of boson creation operators. The critical exponents can be found in terms 
of bulk susceptibilities which can be calculated by conventional Bethe ansatz 
methods. This Luttinger liquid theory was a forerunner of, and is in fact a 
special case of, modern conformal techniques. 

As is discussed further in the next section, the 1-1-1 dimensional Bethe ansatz 
systems are intimately related to two dimensional statistical systems. Essen- 
tially, they are analytic continuations of each other from time to imaginary time. 
Thus techniques developed for one can be readily applied to the other. 

Macroscopic thermodynamic systems at the critical point have been long 
known to be scale invariant. Polyakov|19| first pointed out that they are also 
conformally invariant, that is, invariant under any transformation that is locally 



a scale change plus a rotation. In two dimensions, this constraint dramatically 
reduces possible varieties of critical behavior. The conformal transformations of 
the plane, under which the system must be invariant, correspond to the set of all 
analytic functions, and the generators of these transformations are the elements 
of the so-called Virasoro algebra. Representations of this algebra, correspond- 
ing to possible critical systems, are labelled by a central charge c and other 
parameters usually denoted A^. For the statistical systems, the central charge 
measures the finite-size dependence of the free energy (or the low temperature 
specific heat for the spin chain), and the other parameters measure correlation 
exponents for both systems. It was shown by Friedan, Qiu and Shenker [po| that 
unitarity considerations only allow certain values of c and the A^, and these 
values determines possible excitation energies of the system, in particular the 
low-lying states. But these are precisely the energy levels that can be evaluated 
fairly straightforwardly using the Bethe ansatz, so by this backdoor method we 
can use Bethe ansatz results to learn about correlation functions without trying 
to evaluate matrix elements! A great deal of numerical work has been done to 
classify the conformal theories corresponding to continuum limits of integrable 
spin chain systems. pH 

The Luttinger liquid technique (with some conformal ideas added) has been 
applied successfully to the Hubbard model by Frahm and KorepinQ. This is 
a more complicated situation, in that there are two sets of low energy excita- 
tions, charge and spin density waves, having different speeds. This means that 
conformal methods need to be extended. The basic tool used is the dressed 
charge matrix, a generalized susceptibility, the elements of which can be cal- 
culated using the Bethe ansatz. This immediately gives asymptotic correlation 
functions in terms of coupling strength and magnetic field. Recently, Kane and 
Fisher pOl have applied Luttinger liquid techniques to predicting conductivity 
in a one-dimensional system with barriers. 

4 Transfer matrices and analyticity 

Our growing understanding of the XX Z (and fully anisotropic XYZ) spin one 
half chains has come about in large part through their connections with the 
six (and eight) vertex models. These are classical two dimensional statistical 
models defined on a square lattice, such that on each bond there is an arrow 
which can point either way. The energy, and therefore the statistical weight, 
of each vertex depends on its arrow configuration. In the six vertex model, 
originally conceived to describe hydrogen-bonded crystals (the arrow direction 
indicates the end of the bond where the hydrogen atom sits) electrostatic energy 
considerations dictate that each vertex has two and only two hydrogen atoms 
nearby, so two arrows point in and two out at each vertex. This gives six possible 
vertices, three arrow-reversed pairs. (The eight vertex model allows all in or all 
out arrows.) Taking arrow reversed pairs to have the same energy, there are 



three distinct statistical weights, usually denoted a, 5, and c. The free energy of 
these lattice systems can be found by transfer matrix techniques — the matrix is 
defined between adjacent rows of N vertical bonds, each matrix element giving 
the combined statistical weight of all allowed vertex configurations on the row of 
horizontal bonds sandwiched between the two adjacent rows of specified vertical 
bonds. This is a 2''^ x 2^ matrix with elements having row and column labels 
corresponding to the 2^ possible spin configurations on each of two successive 
rows of vertical bonds. Denoting this matrix by T, the free energy of the N x N 
system is given by TrT^ , so the matrix multiplication sums over all possible 
vertical bond configurations, with the appropriate Boltzmann weighting factor 
supplied by the matrix. For large N, the trace is dominated by the largest 
eigenvalue of T, so this gives the free energy. Finding correlation functions of 
the system requires knowledge of the next leading eigenvalues. 

Now the transfer matrix T operates on a space of N up or down arrows, 
which is the same as the configuration space of the spin-one-half chain Hamil- 
tonian. This was first appreciated by Lieb[p3[, who solved the ice problem by 
using a Bethe ansatz wavefunction for the leading eigenstate of the transfer 
matrix. However, the relationship between T's and Hxxz's is not one-to- 
one. Taking vertex weights with all arrows reversed to be equal, and factor- 
ing out an overall normalization, we still have a fwo-parameter family of ver- 
tex models, with transfer matrices conveniently labelled T(v,"f), where T{v,"f) 
has the same set of eigenstates (but not of course eigenvalues) as the mem- 
ber Hxxzil) of the one-parameter XX Z spin-chain family. That is to say, if 
we take 7 as fixed, there is a one-parameter family of transfer matrices T{v) 
which all commute with Hxxzil), and in fact they all commute with each 
other, since they have the same set of eigenstates. In the notation of Kliimper 
et al.,p3 the ratio of the three vertex statistical weights a : b : c becomes 
sin(7/2 -t- iv) : sin(7/2 — iv) : sin 7, and it can easily be shown that for the 
particular parameter values v = ±17/2, the transfer matrices reduce to simple 
shift operators (because a or 6 become zero), and their logarithmic derivatives 
at these points are just the XX Z Hamiltonian, 

T(=F(i7/2 - v)) = e±'^-''-f^^^^+°('''' (8) 



It might appear at this point that introducing this family of commuting 
transfer matrices does not obviously look like progress in understanding the 
properties of the XX Z spin chain. In fact, though, the existence of the family 
T{v) makes it evident that the XX Z chain is integrable, in the sense of having 
a sequence of higher conserved quantities — i.e., beyond total momentum and 
energy — one could, for example, find the higher logarithmic derivatives of T{v) 
at •(; = ±17/2. However, the real utility of the T{v) arises from their analytic 
properties. As Baxter makes clear throughout his classic book, the analyticity 
of T{v) as a function of v constrains the properties of these systems, and renders 



them solvable, in a truly amazing way. As an example of the connection between 
the Bethe Ansatz equations and analyticity, we mention Baxter's result that 
each eigenvalue A{v) of the transfer matrix satisfies the equation: 

A{v)q{v) = ^{v - i^/2)q{v + ij) + ^{v + i^/2)q{v - 17) (9) 

with $(?;) = (sinhw)^,and q{v) — YVj=i sinh(t; — Vj). The rapidities Vj satisfy 
the Bethe ansatz equations p{vj) = — 1 where 

^ 'i>{v-i-f/2)qiv + ij) 
^^ ' ^v + i^/2)q{v-ij) ^ ' 

These are of course the same equations as in (7) above, but in the context 
of the properties of the transfer matrix eigenvalue as a function of w, we see 
that the zeros of q{v) coincide with zeros of the right hand side of equation (9), 
ensuring that A(t;) has no poles. Furthermore, for the largest eigenvalue, it is 
known that the i^^are all real, from which it follows that A(t;) is analytic and 
nonzero in a strip of width 7 centered on the real axis. This is a very strong 
constraint, and makes it possible to use Cauchy's theorem to derive a nonlinear 
integral equation for the function l/p{x — i7/2), which is exact even for finite 
systems. Although the equation cannot be completely solved, it can be used to 
find the 1/N corrections to the eigenvalues for large systems. This means we 
can evaluate critical exponents analytically [E4| IT] without directly solving the 
Bethe ansatz equations. 

As mentioned earlier, the standard method for evaluating thermodynamic 
properties of Bethe Ansatz systems leads to formidably difficult calculations 
involving large numbers of string-like excitation densities, even the completeness 
of the set of states is sometimes hard to establish. Fortunately, it has been 
realized quite recently that the free energy can also be evaluated more directly 
using the analyticity properties. One needs to evaluate Tr{e~^^^^^). Writing 
this in discretized fashion as 

Tr fe-TT-Hxxz\ (^11^ 

at first glance it appears to be dominated by the largest eigenvalue of the op- 
erator in the bracket. This is of course incorrect, because as N goes to infinity, 
all the eigenvalues degenerate to unity. The trick to getting around this is 
to use crossing symmetry p5| U&,. A neat presentation was given recently by 
Kliimperp5|. From equation (8) above, 

T {i-f/2 - (3/N) .T {-i-t/2 -K (3/N) = e~ « -^^^^ (12) 

Raising this to the N*^ power, one has a product of transfer matrices for a six- 
vertex model having two different vertex strengths, appearing on alternate rows, 
a horizontally striped lattice. To evaluate this in the large N limit, we switch 



to the equivalent product of column transfer matrices. The column transfer 
matrices are all identical, but down each column the vertex strengths alternate, 
as it crosses all the stripes. The eigenvalues of this inhomogeneous transfer 
matrix do not all collapse to unity as N becomes large — the limit is well defined, 
and leads to an integral equation from which one can find the leading eigenvalue. 
Therefore, by this method we can find the thermodynamic behavior of the spin 
chain in a very direct way, without the headache of summing over many different 
string excitations, and concerns about completeness. The new approach proves 
to be very effective computationally. 

5 Inverse-square type systems 

At this point, we back up to our N billiard balls on a line model (meaning 
point bosons), but this time add an inverse-square repulsive potential between 
pairs of billiard balls. This is a famous classical integrable system, solved by 
CalogeroPq, and others, using a Lax pair of matrices. One result that emerged 
was that if initially the balls are far apart with momenta Pi,P2, ■■■,Pn then the 
final state has them far apart with the same set of momenta. This suggested to 
Sutherland p^ the idea of an asymptotic Bethe ansatz for the quantum system. 
The wavefunction is certainly not of the standard Bethe ansatz form of a sum 
of products of plane waves, because the smoothly varying potential energy be- 
tween particles will give functions of smoothly varying wavelength, destroying 
the plane wave property. Nevertheless, the dilute system will be close to plane 
wave form in most of configuration space, and if the scattering phases are fac- 
torizable so that diffraction does not occur, the Bethe ansatz might make sense 
in this dilute limit. Actually, the phase shifting is particular simple for the 1/r^ 
potential, since there is no scale, it is independent of relative momentum, except 
for the sign. To be precise, for a potential of strength 2A(A — l)/r^, the phase 
shift is wXk/ \k\ . This makes the Bethe ansatz equations like (4) above extremely 
easy to solve! Of course, by going to equations like (4), we have switched to 
periodic boundary conditions, but it turns out this can be accomplished by re- 
placing the inverse-square potential by the periodic version, 2A(A — l)/sin^r, 
without spoiling the integrability. In fact, the ground state wavefunction for 
this potential is known exactly. It is rii<7 |sin(a;i — Xj)\ , quite unlike a Bethe 
ansatz wavefunction. Yet, remarkably, the asymptotic Bethe ansatz equations 
give the exact ground state energy of the system even in this manifestly nondi- 
lute limit. They also give all the excitation energies correctly — in other words, 
the full spectrum and the thermodynamics. Note that for 2A = 1,2 or 4 the 
ground state correlation function is exactly the distribution of eigenstates of 
random matrices found by Dysonp7| for orthogonal, unitary and symplectic 
ensembles. 

The high degeneracy of the spectra indicates the presence of new types of 
symmetries psl. There has been much recent work on the algebraic structures 



of these systems, see for example the review of Bernard p9J. An illuminating 
interpretation has been suggested by HaldaneQ, who argues that the magnon 
excitations of the inverse-square spin one-half chain are actually a gas of free 
anyons, specifically spin one-half semions, so the dynamics of the system follows 
entirely from the statistical "interaction" . 

As a simple example of the power of the algebraic approach, we present here 
a new exchange operator formalism invented by Polychronakos |p l| , which gives 
an easy way to prove the integrability of some of these systems. For N particles 
on a line with positions and momenta Xi , pi he introduces an operator Mij that 
permutes the positions and momenta of particles «, j, then defines a generalized 
momentum operator 

TTi = K + « E V{Xi - Xj)Mij (13) 

with V as yet unspecified. He then requires the Hamiltonian to have free-particle 
form in the 7r,'s, 
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gives 
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1>:p^ 




with 
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H-\Y.^1 (14) 

' Pi)M^, + V^^M,, + V^^ - i E V.,kM,,u (15) 



y^,V,k+VJkVk^+Vk^V,, (16) 

and Myfegenerates cyclic permutations, Mijk = MijMjk- 

Requiring the Hamiltonian to have only two-body potentials, V{x) must 
satisfy 

V{x)V{y) + V{y)V{z) + V{z)V{x) = W{x) + W{y) + W{z),x + y + z^O (17) 

The commutator of the tt^'s is 

[^^^^3] = H V,jk{M,jk - Mj^fe) (18) 

At this point, it is very easy to establish the integrability of the Calogero system, 
in other words, to find the N independent constants of motion of the system. 
Assume V{x) = l/x. Then the Hamiltonian has the form 

H=\y^puy.f^^ (19) 

and since on the boson system Mij — 0, this is just the Calogero system. Now 
for V{x) = Ijx, from (16) Vijfe = 0, so from (18) the tt^ all commute, and, 
trivially, so do the symmetric quantities 



/n=E^" (20) 



Assuming now V{x) — I cot ax gives the periodic system. The tt^ no longer all 
commute, but by adding I X]i=^z ^ij to each we generate a new set of /„ that 
do commute, so again we establish integrability with a minimum of effort! 

In fact, this route leads directly to the integrability of the Haldane-Shastry 
spin chain system. All we have to do is to give the bosons a spin of one-half, 
and take the classical (infinitely heavy boson) limit. (This is a mathematical 
model, so we do not worry about spin and statistics.) On these bosons, the 
operator Mij above is essentially equivalent to a spin exchange operator. In the 
classical limit, the bosons form a lattice ground state, and the spin Hamiltonian 
is identical to the Haldane-Shastry model — the momentum term disappears in 
the expression for H above, establishing the integrability of the spin chain [p^. 
The analogous argument applied to the Calogero system of bosons confined by 
an external harmonic potential gives a Haldane-Shastry type system in which 
the spins are located at points corresponding to zeros of an Hermite polynomial. 
Again, the eigenstates are highly degenerate, p^ 

Another inverse-square integrable system that has been extensively investi- 
gated recentlyQ is the Hubbard model with long-range hopping: 

J^ = / ^ tm,nCjna'^nij + U y ^ 'T-mt^mi (21) 

where 

_ it(-l)™-" 

tm.n — T : - (22) 

— sm^ 

This system also has highly degenerate energy levels, given by Bethe ansatz like 
equations which are much more tractable than those for the standard Hubbard 
model, so this is an interesting easy model for exploring properties of one- 
dimensional conductors. 

Finally, wc should at least mention the intriguing connection between inte- 
grability and random matrices, although this will be discussed in more detail 
by Altshuler at this Conference. As stated above, the ground state wavefunc- 
tion for the periodic inverse-square boson system with interaction parameter A, 
Wi^j |sin(a;i — Xj)\ , is for 2A = 1, 2 or 4, the distribution function of eigenvalues 
found by Dyson for orthogonal, unitary and symplectic random matrix ensem- 
bles. Can any meaning be given to the time-dependent correlation functions of 
the boson system in terms of random systems? Simons, Lee and Altshuler [[l2| 
explored how the eigenvalues of a random system (electronic energy levels in 
a disordered conductor) respond to a steadily growing external potential (or 
a twisting boundary condition). For the disordered conductor, they were able 
to use supersymmetry methods developed by Efetov[^ to evaluate eigenvalue 
correlations. They argue that for a potential 

H = Ha + XV (23) 

where Hq is a random potential, containing many nonintcracting electrons, writ- 
ing A^ = 2it, the movement of the eigenvalues as A is increased exactly mirrors 
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the dynamics of the bosons with the inverse-square type interaction. Perhaps it 
is a suitable note to end on that one of the most ordered systems imaginable, an 
integrable system with many higher symmetries reflected in its simple spectrum, 
is intimately connected mathematically to a totally disordered system. 
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